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We theoretically study the electronic properties of BC2N nanoribbons with zigzag edges using 
a tight binding model. We show that the zigzag BC2N nanoribbons have the flat bands and edge 
states when atoms are arranged as B-C-N-C along the zigzag lines. The length of the flat bands in 
the wavevector space depends on the atomic arrangement. This property can be explained by the 
deviation of the linear dispersion of the BC2N sheet from K point of the honeycomb lattice. The 
charge distributions in the edge states depend on the atomic arrangement. We also show that the 
charge distribution of the edge states in zigzag BC2N nanoribbons where the outermost sites are 
occupied with B and N atoms is different from those in conventional graphene zigzag edge. Such 
charge distribution causes different magnetic structures. We investigate the magnetic structure of 
BC2N nanoribbons with zigzag edges using the Hubbard model within a mean field approximation. 
At the zigzag edge where the outermost sites are occupied with B and N atoms, ferromagnetic 
structure appears when the site energies are larger than the on-site Coulomb interaction. 

PACS numbers: 73.20.At, 73.21.Cd, 73.22.-f 



I. INTRODUCTION 

Graphene is an atomically thin carbon sheet in which 
carbon atoms are arranged in a honeycomb lattice. Due 
to its outstanding electronic structure and electron trans- 
port properties, graphene attracts much interests for fu- 
ture electronic devices. Nanocarbon materials which are 
fractions of graphene exhibit different electronic prop- 
erties compared with graphene. Graphene nanoribbons 
are finite width graphene sheets, which are the one of the 
famous examples of nano-carbon materials!^ Recently, 
graphene nanoribbons have been fabricated by e-beam 
lithography 3 , unzipping of carbon nanotubes^* and us- 
ing bottom-up processes^. 

The electronic properties of graph ene n anoribbons 
strongly depend on the edge structures! 1 * 2 * 7 * ^ Graphene 
nanoribbons with armchair edges have the band gaps. 9 
Graphene nanoribbons with zigz ag edg es have the so- 
called flat bands at the Fermi levelPJ^lZEl The states corre- 
sponding the flat bands are loc alized at the zigzag edges. 
They are the edge states^EEl j n the honeycomb lattice, 
there are two inequivalent sites, A- and B-sublattices. 
For the formation of edge states, these sublattice struc- 
tures play decisive roles. The distribution of electronic 
charge of the edge states becomes finite only one sub- 
lattice sites including the outermost sublattice Quite 
recently, the edge states in graphene nanoribbons have 
been confirmed by STM/STS measurement.^ 

On the other hand, hexagonal boron-nitride (BN) 
sheet also shows the honeycomb lattice structure in which 
B and N atoms are arranged in A- and B-sublattice 
sites alternately. Due to their chemical difference be- 
tween B and N, BN sheet has a wide band gap! 11 * 12 ! 
Since B (N) atoms behave as acceptors (donors), boron- 
carbon-nitride-graphene sheet doped with B and N— 
should show interesting controllable electronic proper- 



ties by doping. Ci et al. recently have reported the 
synthesis of hybridized BN and graphene sheets using 
thermal catalytic chemical vapor deposition.^ Such hy- 
bridized BN and graphene sheets show different elec- 
tronic tran sport properties compared with graphene and 
BN sheetsFGaHl The electronic properties of BN and 
graphene hybridized nanoribbons have been investigated 
by several authors ff^H^ In such nanoribbons, ferrimag- 
netic ordering along the edges were predicted P2HHI Re- 
cently, Kaneko et al. reported that the appearance of flat 
bands and edge states in zigzag BCN nanoribbons where 
the outermost C atoms of graphene nanoribbons are re- 
placed with B and N atoms alternately.^ They reported 
that the distributions of charge and spin densities in the 
edge states are different from those in the edge states of 
zigzag graphene nanoribbons, i.e., the electronic charge 
distributes over both sublattices and the spin density can 
become ferromagnetic at the edges. 

BC2N sheet is organic analogous of graphene, which 
can be regarded as one of example of boron-carbon- 
nitridc. Graphite-like BC2N has been synthesized using 
chemical vapor depositions of boron trichloride, BCI3, 
and acetronitrile, CH 3 CN as the starting materials! 20 * 21 ! 
The electronic properties of BC2N sheets depend on the 
atomic arrangement! 2 ^" 2 ^ The electronic properties of 
nanoribbons made with BC2N have been also investi- 
gated by several authorsP^U l u e f a i reported that 
BC2N nanoribbons can have the magnetization even for 
armchair edge nanoribbons J2512S1 Xu et al. showed that 
BC2N nanoribbons with zigzag edges have the linear dis- 
persion when the atoms are arranged C-B-N-C in the 
transverse direction However, there are no reports on 
the presence of the flat bands and edge states in BC2N 
nanoribbons. The purpose of this paper is to explore the 
flat bands and edge states in BC2N nanoribbons. 

In this paper, we will investigate the electronic proper- 
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FIG. 1. Models of BC2N nanoribbons handled in this study. In this figure, the black and white circles represent B and N 
atoms, respectively. C atoms are located at the empty vertices. The dotted rectangles represent the unit cells and the shaded 
regions represent the unit cell of BC2N sheets. 



ties of BC2N nanoribbons with zigzag edges using a tight 
binding model. We will show that zigzag BC2N nanorib- 
bons have the flat bands and edge states when atoms are 
arranged as B-C-N-C along the zigzag lines. The charge 
distribution of the edge states in zigzag BC2N nanorib- 
bons is different from those in conventional graphene 
zigzag edge. Corresponding magnetic structures are in- 
vestigated using the Hubbard model within a mean field 
approximation. 

This paper is organized as follows: In Sec. II, the mod- 
els of BC2N nanoribbons and the Hamiltonian of the sys- 
tem are presented. The spin independent results within 
the tight binding model and the spin dependent results 
within the Hubbard models are presented in Sec. Ill A 
and III B, respectively. The discussion and short sum- 
mary of the paper are given in Sec. IV and V, respectively. 
In Appendix A, the results with B-B and N-N bondings 
are summarized. 



II. MODELS AND METHODS 

We shall consider six different structures of BC2N 
nanoribbons with zigzag edges as shown in Fig. [l] In 
this figure, B and N atoms, respectively, are indicated 
by the black and white circles, and C atoms are lo- 
cated the empty vertexes. It should be noted that atoms 
are arranged as B-C-N-C along the zigzag line in these 
BC2N nanoribbons. The dotted rectangles represent the 
unit cells and the shaded regions represent the unit cell 
of BC2N nanoribbons and the shaded regions represent 
BC2N sheets which consist of the nanoribbons. Let N 
be a number of the zigzag lines. In this study, we shall 
restrict ourself to consider even N cases for simplicity. 
Because these bondings are not energetically favorable, 
we do not consider the BC2N nanoribbons with B-B or 
N-N bondings. The results of such nanoribbons are sum- 
marized in Appendix A and the discussion about this 
issue is given in Sec. II C. 



The Hamiltonian of the system within the tight bind- 
ing model of 7r-electrons is given by 

n = J2 EiC i c i ~ H hAcj' (!) 
» (id) 

where Ei is an energy of tt electron at the site i, c] and Cj 
are the creation and annihilation operators of electrons 
at the lattice site i, respectively, stands for sum- 

mation over the adjacent atoms and tij is the hopping 
integral of n electrons from jth atom to ith atom. Ei 
are the site energies, Eb, Eq and E^, at the B, C and N 
sites, respectively. Following to the Yoshioka et aL, we 
shall assume that the hopping integrals are constant re- 
gardless of the atoms, tij = t, E^ = —Eb and Ec = Op^ 
In Ref. [21], E-g/t has been estimates as 0.8 ~ 0.9. How- 
ever, the hopping integral, t = 3.0 eV, used in Ref. [H] 
seems to be overestimated compared with that deter- 
mined by the first principles calculations within the local 
density approximation (LDA), £ = 2.6 — 2.7 eVP 1 On 
the other hand, the Shubunikov-de Haas meas urements^ 
and the first-principles GW calculations^*^ g ave much 
larger values such as t = 3.4 ~ 3.5 eV. Furthermore, the 
site energy, Eb, also have ambiguity. In BN sheet, the en- 
ergy gap is given by 2Eb- The band gap was obtained 3.9 
eV and 5.4 eV by the first-principles calculations within 
LDA and GW approximation.^ Ribeiro and Peres es- 
timated as i = 2.33 eV and Eb — 1-96 eV by fitting 
the to the band structure of BN sheet within the tight 
binding model to that of the first-principles calculation 
within the generalized gradient approximation (GGA) 1^3 
Zheng et al. proposed that t = 2.621 eV, E B = 1.894 
eV and the overlap integral s — 0.0154 eV by fitting to 
the band structure of BN nanoribbon within the tight 
binding model to that of the first-principles calculation 
within LDA.^But, Eb can vary about 2 eV depending 
on their environment.^ To investigate the dependence of 
Eb on the electronic properties, therefore, we shall con- 
sider E B /t = 0.7, 1.0 and 1.3. 

In order to discuss the magnetism in BC2N nanorib- 
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FIG. 2. The band structures of BC2N nanoribbons for TV = 10 and for the model-1 (a), model-2 (b), model-3 (c), model-4 
(d), model-5 (e) and model-6 (f). The results of Es/t — 0.7 (i), 1.0 (ii) and 1.3 (iii) are presented in the left, middle and right 
panels, respectively. The downward arrows indicate the position of k — 2n/3. 



bons, we shall use the Hubbard model. The Hamiltonian 
of the Hubbard model can be obtained by including the 
on-site Coulomb interaction to Eq. ([I]): 

H = ^2 Eicj ^a - ^ U,A ia Cj,cr + ^ U i n i,t n i,h ( 2 ) 

where U% is the on-site Coulomb interaction at «th site 
and rii jC = c\ a Ci tCr (a =t, I) is the number operator of 
electron at it site with spin a . In this study, we shall 
assume that Ut is also constant, [/PS We adopted the 
mean field approximation for this Hamiltonian in order 
to solve the problem. The spin expectation values at ith 
sites are defined as Si tZ — ((n^f) — (7i^))/2 

The rang e of m agnitude of U is characterized by that 
of graphene! 30 * 36 * In irans-polyacetylene which can be re- 
garded as one-dimensional limit of graphene na noribb ons 
with zigzag edges, U was estimated as 3.0 eV! 37 | 3S | Ac- 
cording to the first-principles calculations within the lo- 
cal spin density approximation (LSDA) and G GA, U/t 
was estimated as 0.9 and 1.3, respectively ! 39 * 40 ! The mag- 
netic structures in BCN nanoribbons within LSDA can 



be well reproduced in the Hubbard model.^ In the pre- 
vious study, < Eq ft < 2 and < U/t < 2 were used 
and the results within LSDA can be well interpreted. In 
this paper, we shall consider < U/t < 2. 



III. RESULTS AND DISCUSSION 
A. Spin independent results 

In Fig. [2j the dependence of E-q on band structures of 
BC2N nanoribbons for = 10 and for the model-1 (a), 
the model-2 (b), the model-3 (c), the model-4 (d), the 
model-5 (e) and the model-6 (f) are shown. The results 
of Es/t = 0.7, 1.0 and 1.3 are presented in the left (i), 
middle (ii) and right panels (iii) , respectively. The down- 
ward arrows indicate the position of k = 2ir/3, i.e., the 
projection of the K point of two dimensional graphene 
to the one dimensional system. For E-g/t = 0.7 cases, 
we have found the flat bands at E = in the all BC2N 
nanoribbons. The length of the flat bands in wavevector 
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FIG. 3. Same figure 

space seems depending on the atomic arrangements. For 
Es/t = 1.0 cases shown in the middle panels, we have 
observed the flat bands at E = but the length of the flat 
bands have been changed in some models, too. With in- 
creasing i?B, the length of the flat bands increases for the 
model-2 and -3 nanoribbons, while the length decreases 
for the model-5 and -6 nanoribbons. On the other hand, 
for the model-1 and -4 nanoribbons, the length of the 
flat bands is not sensitive to the magnitude of E#. For 
Es/t = 1.3 cases, these changes in band structures can 
be clearly seen. For the model-2 and -3 nanoribbons, the 
length of the flat bands increases but the degeneracies of 
flat bands are lifted. On the other hand, for the model-5 
and -6 nanoribbons, the flat bands become much shorter. 
These results might indicate that there are threshold of 
Eb to have the flat bands. 

In order to see the width dependence of the band struc- 
tures, we shall consider the much wider BC2N nanorib- 
bons. In Fig. [3j the similar figures as Fig. [2] for N = 30 
are present. We have found that there are flat bands 
at E = for all nanoribbons independent of Eb- Fur- 
thermore, we can clearly see that the length of the flat 
bands increases for the model-2 and -3 nanoribbons but 
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as Fig. 2 for N = 30. 

decreases for the model-5 and -6 nanoribbons with in- 
creasing Eb- The length of the flat bands in the model- 1 
and -4 are independent of Eb- Therefore, the length of 
the flat bands depends on the the atomic arrangement 
and Eb, and there is not a threshold value for the pres- 
ence of the flat bands. It should be emphasized that we 
have confirmed that the flat bands are absent if atoms are 
not arranged as B-C-N-C along the zigzag lines. There- 
fore, we can conclude that B-C-N-C arrangement along 
the zigzag line is necessary to obtain the flat bands. 

In order to understand the difference in the length of 
the flat bands in the wavevector space, we shall con- 
sider the band structure of BC2N sheet. In Figs. [4] (a), 
schematic illustrations of BC2N sheets are present. The 
shaded region indicates the unit cell for the band calcu- 
lations. It is noted that the model-1 nanoribbon consists 
of BC2N sheet shown in Fig. [4] (a)-(i) while the model- 
2 and -3 nanoribbons consist of BC2N sheet shown in 
Fig. [4] (a)-(ii). Figure [4] (b)-(i) and (ii) show the band 
structures of BC2N sheet for £?bA = 1-0 shown in Fig. [4] 
(a)-(i) and (ii), respectively. In Fig. [4] (b)-(iii), the band 
structure of graphene in 2 x 2 supercell is presented as 
a reference. For BC2N sheet shown in Fig. [4] (a)-(i), the 



FIG. 4. (a) Schematic illustration of BC2N sheets presence 
of hexagonal symmetry. BC2N sheet shown in panel (i) con- 
sists the model-1 BC2N nanoribbon and BC2N sheet shown 
in panel (ii) consists the model-2 and 3 BC2N nanoribbons. 
(b) The band structures of BC 2 N sheet for Es/t = 1.0. Panel 
(i) and (ii) show the band structures of BC2N sheets shown 
in panel (i) and (ii) of Fig. (a), respectively. In panel (iii) 
corresponding band structure of graphene is shown, (c) The 
Brillouin zone of the honeycomb lattice (the dashed lines) and 
its projection to the one dimensional Brillouin zone of zigzag 
nanoribbons. 



linear dispersion is realized at the K point. This fea- 
ture is quite similar to that of graphene but the group 
velocity around the Fermi level is slightly smaller than 
that of graphene. On the other hand, for BC2N sheet of 
Fig. [3] (a)-(ii), the position of linear dispersion deviates 
from the K point toward the T point. The deviation of 
the linear dispersion increases with increasing Eb- This 
result is quite similar to the result calculated by Liu et 
al. using the first principles calculations.^ As discussed 
below, the position of the Dirac points plays decisive role 
to determine the length of flat bands. 

Figure [4] (c) shows the Brillouin zone of the honey- 
comb lattice and its projection to the one dimensional 
Brillouin zone of zigzag nanoribbons. The projection of 
the K point of the honeycomb lattice, i.e., the Dirac point 
of graphene and BC2N sheet shown in Fig. [I] (a)- (i) is lo- 
cated at k = 2n/3. Then, the model-1 nanoribbons have 
the flat bands in < k < 2ir/3 which is same region 
as graphene nanoribbons. When the Dirac points devi- 
ate from the K point toward the T point as indicated 
by the arrows, the projection of the Dirac points shift 
from k = 27r/3 toward k — tt as shown in the bottom of 
Fig. [4] (c). So, the length of the flat bands increases. As 
the result, the difference in the band structures of BC2N 
sheets might lead to the difference in the length of the 
flat bands in wavevector space. 



FIG. 5. (a) Schematic illustration of BC2N sheets absence 
of hexagonal symmetry. BC2N sheet shown in panel (i) con- 
sists the model-4 BC2N nanoribbon and BC2N sheet shown 
in panel (ii) consists the model-5 and 6 BC2N nanoribbons. 
(b) The band structures of BC 2 N sheet for E B /t = 1.0. Panel 
(i) and (ii) show the band structures of BC2N sheets shown 
in panel (i) and (ii) of Fig. (a), respectively. In panel (iii) 
corresponding band structure of graphene is shown, (c) The 
Brillouin zone of the rectangular lattice. 



Next, we shall consider the BC2N sheet shown in Fig. 
[5] The model-4 nanoribbon consists of BC2N sheet 
are presented in Fig. [5] (a)-(i) while the model-5 and - 
6 nanoribbons consist of BC2N sheet are shown in Fig. 
[5] (a)-(ii). In Fig. [5] (a), the band structures of BC 2 N 
sheets for E^/t = 1.0 are shown in Fig. [2] (b). The band 
structure of graphene in the same unit cell is shown in 
Fig. [5] (b)-(iii). The corresponding Brillouin zone is in- 
dicated in Fig. [5] (c) . The vertical dotted line in Fig. [2] 
(b) indicates the position of the K point of honeycomb 
lattice. For BC2N sheet shown in Fig.[5](a)-(i), the linear 
dispersion is located at the same point of the graphene. 
As displayed in the lower part of Fig. [5] (c), the pro- 
jection of the K point of honeycomb lattice is located 
at k = 27r/3 in the one dimensional Brillouin zone of 
nanoribbons. Therefore, the length of the flat bands in 
the wavevector space is same as those in the graphene 
nanoribbons. On the other hand, the point of linear dis- 
persion shifts toward the T point in BC2N sheet shown 
111 Fig. [5](a)-(ii). In this case, the position of the linear 
dispersion in the one dimensional Brillouin zone also shift 
toward k = 0. Then, the length of the flat bands of the 
model-5 and -6 BC2N nanoribbons is shorter than those 
of graphene. Furthermore, the length of flat bands of 
other BC2N nanoribbons also vary in a similar manner 
as discussed Appendix A. To conclude, we have found 
that the length of the flat bands in BC2N nanoribbons 
is determined by the position of the Dirac points in the 
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FIG. 6. Calculated local density of states (LDOS) for Es/t = 1.0 and N = 10 with schematic illustrations for the model-1 
(a), model-2 (b), model-3 (c), model-4 (d), model-5 (e) and model-6 (f). As the reference, LDOS of graphene nanoribbons in 
doubled unit cell is presented in (g). In the left side of each panel, the corresponding schematic illustration is presented. The 
magnitude of LDOS is indicated by the area of the circles. 



Brillouin zone. 

Next, we move to the charge distribution in zigzag 
BC2N nanoribbons. In Fig. |6j the local density of states 
(LDOS) at E = for several structures with E^/t = 1.0 
are shown by the circles. In this figure, the area of the 
circles is proportional to the magnitude of the LDOS at 
each site. In the left side of each panel, the corresponding 
schematic illustration is presented. As the reference, the 
LDOS at E = graphene nanoribbon is shown in Fig. [6] 
(g). In the model-1 and -4 nanoribbons shown in Figs.p^ 
(a) and (d), the upper and lower edge are different each 
other. In the following, we shall call the upper (lower) 
edge as B, N side (C side) edge. The electronic charge 
is localized at the BC2N nanoribbons edges, showing the 
presence of the edge states. As discussed below, the edge 
states in BC2N nanoribbons are different from those in 
conventional graphene nanoribbons. 

In the model-1, the charge distributions along the B, 
N side and C side edges are different each other. The 

model-2 



k = n k=0.8n fc=0.67i /(=0.4ji 




FIG. 7. The charge distribution of the edge states with 
E > in the model-2 BC2N nanoribbons for N = 10 and 
Es/t = 1.0. The areas of the circles indicate the magnitude 
of the charge densities. 



charge distribution at the C side edge is similar to that 
at the conventional zigzag edge of graphene nanorib- 
bons, while the charge of the edge states at the B, N 
side edge distributes the both sublattice sites. Recently, 
Kaneko et al. have shown that the edge states in zigzag 
graphene nanoribbons are robust on the substitution of 
outermost C atoms with B and N atoms alternately.^ 
However, such substitution causes change in charge dis- 
tribution, i.e., the sublattice structure is brokenPSJ The 
edge states at the B, N side edge is similar to those discov- 
ered by Kaneko et c/iP^ In the model-2 nanoribbon, the 
charge distribution of the edge states is similar to that of 
graphene nanoribbons but the sublattice structure is bro- 
ken inside the nanoribbons. In the model-3 nanoribbon, 
the charge distributes over both sublattice sites, show- 
ing the similarity of those discovered by Kaneko et aP^ 
The charge distribution in the model-4 BC2N nanoribbon 
is quite similar to that in the model-1 BC2nanoribbon. 
The charge distributions in the model-5 and -6 BC2N 
nanoribbons are similar to those in the model-2 and - 
3 BC2N nanoribbons, respectively. However, the LDOS 
inside of nanoribbons in the model-5 and -6 are smaller 
than those in the model-2 and -3 nanoribbons. This dif- 
ference is caused by the difference in the length of the 
flat bands in the wavevector space. 

In the model-2 nanoribbons, there are other edge states 
even E ^ 0. As shown in Fig. [2] (b) and [3] (b), we have 
found width independent energy bands around E/t ~ ±1 
in k/w ~ 0.2 — 1. The wavevector dependence of the 
charge distribution for N = 10 E-g/t = 1.0 of the edge 
states with E > is shown in Fig. [7] In this figure, 
the areas of the circles indicate the magnitude of the 
charge densities. We can clearly see the formation of the 
edge states even E 7^ 0. At k — tt, the electrons are 
mainly accumulated in the outermost C atoms and the 
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FIG. 8. The spin dependent band structures of BC2N nanoribbons for N = 10 and E^/t — 1.0, and for the model-1 (a), 
model-2 (b), model-3 (c), model-4 (d), model-5 (e), model-6 (f). The results of U/t = 0.5 (i), 1.0 (ii) and 1.5 (iii) are presented 
in the left, middle and right panels, respectively. 



next outermost B atoms, while the electronic charge is 
absent at the next outermost N atoms. With decreas- 
ing of k, the electronic charge at the next outermost N 
atoms increases, resulting in the decrease in the energy. 
For E < states, on the other hand, the electrons are 
mainly accumulated in the outermost C atoms and the 
next outermost N atoms. 



B. Spin dependent results 

Next, we shall consider the magnetism in BC2N 
nanoribbons. Figure [8] shows calculated band structures 
of BC 2 N nanoribbons for E B /t = 1.0 and N = 10. The 
results of U/t = 0.5, 1.0 and 1.5 are presented in the 
left, middle and right panels, respectively. In Figs. [8] (a) 
and (d), the solid and dashed lines represent the energy 
bands of up and down spins, respectively. It should be 
noted that the band structures of up and down spins of 



the model-2, -3, -5 and -6 nanoribbons shown in Fig. [8] 
(b), (c), (d) and (c) are double degenerate. 

We found that the degeneracy of flat bands are lifted 
due to the presence of the Coulomb interaction indepen- 
dent of the atomic arrangement. In the model-2 nanorib- 
bons, the degenerate flat bands of E ^ are also lifted. 
These BC2N nanoribbons are semiconductor with direct 
gap. The band gap increases with increasing U/t. The 
position of valence tops and conduction bottoms are lo- 
cated around k = OAir ~ 0.77T depending on the atomic 
arrangement, showing good agreement with the position 
of projection of the Dirac point. 

Corresponding spin density distributions for U/t = 0.5 
(middle) and 1.5 (right side) are present in Fig. [9] with 
their schematic illustrations (left side). In this figure, the 
up (down) spin densities are denoted by the solid (open) 
circles whose areas are proportional to the magnitudes 
of the spin densities. In the model-1 nanoribbon, the 
ferrimagnetic order is observed at the C terminated edge 
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FIG. 9. Calculated spin densities for E&/t = 1.0 and N = 10 with schematic illustrations (left side), and for the model-1 
(a), model-2 (b), model-3 (c), model-4 (d), model-5 (e) and model-6 (f). The results for U/t — 0.5 and 1.5 are presented in 
the middle and right side in each panel. The up (down) spin densities are denoted by the solid (open) circles whole areas are 
proportional to the magnitudes of the spin densities. 



while the ferromagnetic order is observed in B and N ter- 
minated edge for U/t — 0.5. For U/t = 1.5, on the other 
hand, the ferromagnetic orders are observed at the both 
edges. In the model-2 nanoribbon, the ferrimagnetic or- 
ders are observed at the both edges independent of the 
magnitude of U. However, the magnetic structures in the 
model-3 depend on the magnitude of U while the model-2 
and -3 nanoribbons are consisted with same BC2N sheet. 
We observed ferromagnetic order at the both edges for 
U/t = 0.5 but it turns into ferrimagnetic with increas- 
ing in U. Therefore, we found that the magnetic struc- 
tures in BC2N nanoribbons depend on both the atomic 
arrangements and the magnitude of U. 

It should be noted that the ferromagnetic orders were 
observed when the outermost site are occupied by B and 
N atoms alternately and U/t = 0.5. At these edges, 
LDOS shown in Fig. [6] distributed over both sublattice 
sites. For the ferrimagnetic orders even for U/t = 0.5, 
on the other hand, LDOS at the next outermost sites 
are absent. Therefore, the charge distributions seem to 
determine the magnetic structure in BC2N nanoribbons. 
These behaviors are quite resemble to the previously re- 
ported magnetic orders in BCN nanoribbons where the 
outermost C atoms of graphene nanoribbons are replaced 
by B and N atoms alternately. 

To confirm this conjecture, we shall see the spin distri- 
butions in the model-4, 5 and 6 nanoribbons. The mag- 
netic structure of the model-4 nanoribbons for U/t = 0.5 
is the ferrimagnetic at C side edge but ferromagnetic at 
B,N side edge. For U/t = 1.5, the magnetic structures 
at both edges are the ferrimagnetic. The magnetic struc- 
tures of the model-4 nanoribbon shown in Fig. [9] (d) are 
quite resemble to those of the model-1, showing good 
agreement of the conjecture. In the model-5 nanorib- 
bons, we observed the ferrimagnetic structures indepen- 



dent of the magnitude of U. In this nanoribbons, the out- 
ermost sites are occupied by C atoms and LDOS at the 
next outermost sites are absent. In the model-6 nanorib- 
bons, on the other hand, the ferromagnetic orders are 
observed for U/t — 0.5 but they turn into ferrimagnetic 
with increasing U. The outermost sites of the model-6 
nanoribbons are occupied with B and N atoms and the 
electronic charges distribute over both sublattice sites. 
We found that the magnetic structures in the model-5 
and -6 nanoribbons satisfy the conjecture. 

We shall consider Eb and U dependence of magnetic 
structures in order to to clarify when change in magnetic 
structure takes place. The dependence of U of the spin 
densities are summarized in Fig. |10| In this figure, spin 
expectation values at the outermost sites and the next 
outermost sites are presented in the top and bottom pan- 
els, respectively. The results of E^/t — 0.7, 1.0 and 1.3 
are represented by the solid, dotted and dashed lines, re- 
spectively. It should be emphasized that obtained spin 
expectation values in B and N atom are same as we have 
seen in Fig. [9] Therefore, we plotted the spin expectation 
values in C atoms and B atoms. 

At the B, N side of model-1 nanoribbons, we obtained 
positive s z independent of Eb at the outermost and next 
outermost sites for sufficiently small U. With increasing 
U , s z at the outermost site monotonically increases while 
that at the next outermost site monotonically decrease, 
showing that the magnetic structures change from the 
ferromagnetic into ferrimagnetic order depending on Eb . 
At the C-side of model-1, on the other hand, s z of the out- 
ermost and next outermost sites are negative and zero, 
respectively, for sufficiently small U . With increasing in 
U, s z of the outermost monotonically decreases but that 
of next outermost sites increase, i.e., the ferrimagnetic 
order is observed. Since B and N atoms are arranged in 
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FIG. 10. U dependence of spin expectation values at the outermost sites (i) and the next outermost sites (ii) at B,N side in 
the model-1 (a), C side in the model-1 (b), model-2 (c), model-3 (d), B,N side in the model-4 (e), C side in the model-4 (f), 
model-5 (g) and model-6 (h) BC2N nanoribbons. The solid, dotted and dashed lines represent the results for Es/t = 0.7, 1.0 
and 1.3, respectively. 



the sublattice sites which do not belong to the outermost 
sites, the Eb dependence of s z is small at the C-side 
compared with those at the B, N-side. In the model- 
4 nanoribbons, calculated s z at both edges show quite 
similar dependence on U to those in the model- 1 nanorib- 
bons since the electronic properties of these nanoribbons 
around E = resemble between them. 

In the model-2 nanoribbons, s z are almost independent 
of Eb except U/t < 0.2 regime. As we discussed before, 
the degeneracy of flat bands are lifted with inceasing of 
Eb, resulting in decreasing of strength of the electron 
correlation effect. Therefore, s z becomes smaller with 
increasing Eb- Similar effect can be found in the model- 
3 nanoribbons. Especially, s z vanishes for U/t < 0.1 
in the model-3 nanoribbons with E-a/t = 1.3. With in- 
creasing U, s z at the next outermost site once increases 
but turns to decrease to negative value at certain value 
of U. Then, the ferrimagnetic order is constructed. In 
the model-5 and -6 nanoribbons, obtained s z shows sim- 
ilar dependence on U to those in the model-2 and -3 
nanoribbons, respectively, but the magnitudes of s z be- 
comes smaller since the length of flat bands are shorter 
than those in the model-5 and -6 nanoribbons. As we 



mentioned before, the degeneracies are recovered with in- 
crease of the ribbon width, suggesting that the strength 
of electron correlation effect should increases. It should 
be emphasized that we confirmed that such strong Eb 
dependence for small U will decrease with increasing the 
ribbon width. 



IV. DISCUSSION 

We did not consider BC2N nanoribbons in which atoms 
are arranged as B-C-N-C along zigzag lines with B-B and 
N-N bonds. The results of BC2N nanoribbons with B-B 
and N-N bonds are summarized in Appendix A. We con- 
firmed that such BC2N nanoribbons also have flat bands 
at the Fermi level but the length of the flat bands of such 
BC2N nanoribbons also different from that in conven- 
tional zigzag graphene nanoribbons. It should be noted 
that the change is the length of the flat bands can be 
understood as the change in the Dirac point. The mag- 
netic structures in these BC2N nanoribbons were calcu- 
lated using the Hubbard model within the mean field ap- 
proximation. For the magnetic structures in these BC2N 
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(a) model-7 



(c) model-9 



LDOS 




FIG. 11. The models of BC 2 N 
nanoribbons handled in the appendix 
(left sides) and calculated LDOS at 
E = for E B /t = 1.0 (right sides). 
The dotted rectangles represent the 
unit cells and the shaded regions rep- 
resent the unit cell of BC2N sheets. 



nanoribbons, the atomic arrangement of the outermost 
sites and the ratio of the site every to the on-site Coulomb 
interaction also play decisive role as discussed above. 

Dependence of atomic arrangement on the energetics 
of BC2N sheet has been investigated by Liu et al. and 
Azevedo using the first-principles calculations! 22 * 24 ! Ac- 
cording to their study, BC2N sheets stabilize with in- 
creasing B-N and C-C bondings, suggesting that BN is- 
lands and graphene islands hybridized structure is most 
stable. Such tendency was confirmed in the arbitrary 
stoiciometric BCN, i.e., B 2: C !/ N 2 sheets, by the several 
authorsPH^ Therefore, the BC2N nanoribbons treated 
in this paper are not stable. 

While the atomic structures of synthesized BC2N was 
not identified, Watanabe et al. reported the bonding 
characteristic of BC2N tho films by means of x-ray pho- 
toelectron s pec troscopy measurement of chemical shift of 
Is electrons.^ According to their study, BC2N thin films 
have significant B-C and C-N bondings, i.e., all B, C and 
N atoms bond with one another and mixed atomically.^ 
Therefore, the synthesized BC2N might have the atomic 
arrangement considered in this paper. The synthesized 
graphitic BC2N showed p-doped semiconducting charac- 
ter with band gap about 2 eV by STM and photolumines- 
cence measurement.^ However, BC2N sheets considered 
in this paper are metallic as shown in Figs. [4] (b) and [5] 
(b) . The tight binding model cannot describe the charge 
transfer as discussed in Ref. due to its simplicity. 
The effect of char ge t ransfer might become important in 
BCN nanoribbons PS Therefore, the discussion on the en- 
ergetics and electronic structures of BC2N nanoribbons 
is desired. 

We also performed the first principles calculations 
based on the density functional theories within the pro- 
jector augmented wave ( PAW) methooP and LDAEsI im- 
plemented in vasp code.^ESl The BC 2 N nanoribbons 
shown in Fig. [T] have the flat bands BC2N nanoribbons 
when the outermost atoms are terminated by single H 
atoms. According to the results within the first-principles 
calculations, the on-site energy, E^, seems to be larger 
than the hopping integral, t. However, these issues are 



beyond the scope of the present paper and will be re- 
ported in future publications. 



V. SUMMARY 

We theoretically studied the electronic and magnetic 
properties of BC2N nanoribbons with zigzag edges us- 
ing the tight binding model and the Hubbard model. 
We showed that zigzag BC2N nanoribbons have the flat 
bands and edge states when the atoms are arranged as 
B-C-N-C along the zigzag lines. The length of the flat 
bands depends on the atomic arrangement, which can be 
explained by the shift of the Dirac point in BC2N sheet 
from the K point of the honeycomb lattice. The charge 
distribution also depends on the atomic arrangement. 
When the outermost sites are occupied by C atoms, the 
charge distributions of the edge states in zigzag BC2N 
nanoribbons are resemble to that of conventional edge 
states. When the outermost sites are occupied by B and 
N atoms, on the other hand, the charge of the edge states 
distribute over both sublattice sites, which are different 
from those in conventional graphene zigzag edge. We 
also showed that such charge distribution causes differ- 
ent magnetic structures. When the outermost sites are 
occupied by C atoms, we observed the ferrimagnetic or- 
der at the edges independent of the magnitude of the on- 
site Coulomb interaction, which are resemble to that of 
conventional edge states at the zigzag graphene nanorib- 
bons. At the zigzag edge where the outermost sites are 
occupied with B and N atoms, on the other hand, the 
ferromagnetic structure appears when the site energies 
are larger than the on-site Coulomb interaction and the 
magnetic structures turn into the ferrimagnetic with in- 
creasing the on-site Coulomb interaction. 
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FIG. 12. The band structures of BC2N nanoribbons for TV = 10 and for the model-7 (a), model-8 (b), model-9 (c) and 
model-10 (d). The results of Es/t = 0.7 (i), 1.0 (ii) and 1.3 (iii) are presented in the left, middle and right panels, respectively. 
The downward arrows indicate the position of k = 2tt/3. 
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FIG. 13. Same figure as Fig. p] for N = 30. 
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Appendix A: BC2N nanoribbons with B-B and N-N 
bondings 



In this appendix, we shall present the results of the 
BC2N nanoribbons with B-B and N-N bondings. In the 
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FIG. 14. (a) The band structures of BC 2 N sheet which 
consists the model-7 and 8 BC2N nanoribbons. (b) The band 
structures of BC2N sheet which consists the model-9 and 10 
BC2N nanoribbons. The vertical dotted line represents the 
position of the K point of the honeycomb lattice. 



left sides of Figs. 11 (a-d), the schematic illustrations of 
BC2N nanoribbons with B-B and N-N bondings are pre- 
sented. The calculated band structures for N = 10 and 
30 are summarized in Figs. 12 and 13 respectively. We 



found that the length of the flat bands in the model-7 and 
-8 nanoribbons are shorter than that of zigzag graphene 
nanoribbons, while the flat bands in the model-9 and 
10 nanoribbons are longer than that of zigzag graphene 



nanoribbons. It should be noted that the flat bands of 
the model-9 and 10 nanoribbons for Es/t > 1.0 cover 
the whole Brillouin zone. Calculated LDOS at E = 



are given in right sides of Figs. 11 (a-d) with the magni- 
tude of LDOS are represented by the area of the circles. 
The flat bands correspond to the edge states. We found 
that the electric charge distributes with similar manner 
to those shown in Fig. [6| 

As we discussed above, we shall consider the band 
structures of BC2N sheets. Figure 14 (a) show calcu- 
lated band structures of BC2N sheets which consist the 
model-7 and 8 BC2N nanoribbons. Corresponding Bril- 
louin zone is shown in Fig.|4](c). With increasing Eb, the 
Dirac point deviates from K point toward the M point. 
Therefore, the Dirac point in the one dimensional Bril- 
louin zone shifts toward the T point, resulting in shorter 
flat bands. 

Figure 14 (b) show calculated band structures of BC2N 
sheets which consist of model-9 and 10 BC2N nanorib- 
bons. Corresponding Brillouin zone is shown in [5] (c) . In 
this figure, the vertical dotted line represents the position 
of the K point of the honeycomb lattice. In this case, the 
Dirac point is deviated from the K point of honeycomb 
lattice toward X point for E-g/t — 0.7. Therefore, the 
length of flat band becomes longer than that of zigzag 
graphene nanoribbons. For E^/t = 1.0, the dispersion 
becomes massive at the X point. Then, the flat bands 
cover the whole Brillouin zone. With increasing Eb, we 
observed the Dirac cone but the Dirac point are located 
on the X-K line. In the one dimensional Brillouin zone, 
there are states with E = at k = ir and the flat bands 
remain covering the whole Brillouin zone. Therefore, the 
length of flat bands depend on the atomic arrangement 
as discussed before and the change in the length of flat 
bands can be understood as the shift of the Dirac point 
in BC2N sheet from the K point of honeycomb lattice. 
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